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Abstract 

In the framework of heavy quark effective theory we use QCD sum rules to 
calculate the masses of the cs (0 + ,l + ) and (l + ,2 + ) excited states. The results are 
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consistent with that the states D s j (2317) and D s j (2460) observed by BABAR and 
^ ■ CLEO are the + and 1 + states in the j\ = i + doublet. 

PACS number: 12.39.Hg, 13.25.Hw, 13.25.Ft, 12.38.Lg 

a: 

p_; 1 Introduction 

Recently BaBar Collaboration announced a positive-parity narrow state with a rather 
low mass 2317 ± 3 MeV in the _D+(1969)7r channel [1], which was confirmed by CLEO 
[2] and BELLE [3] later. Because of its low mass and decay angular distribution, its 
J p is believed to be + . In the same experiment CLEO [2] observed a state at 2460 
MeV with the possible spin-parity J p = 1 + in the D*n channel. BaBar [1] also found a 
signal near there. Since these two states lie below DK and D*K threshold respectively, 
the potentially dominant s-wave decay modes _D s j(2317) — > D S K etc are kinematically 
forbidden. Thus the radiative decays and isospin- violating strong decays become favorable 
decay modes. The later decay goes in two steps with the help of virtual D s r] intermediate 
states, _D s j(2317) — > D s r] — > D s tt° where the second step arises from the tiny 77-71"° isospin- 
violating mixing due to m u 7^ ma- 

The experimental discovery of these two states has triggered heated debate on their 
nature recently in literature. The key point is to understand their low masses. The 
-D s j(2317) mass is significantly lower than the values of + mass in the range of 2.4 — 2.6 
GeV calculated in quark models of [4]. The model using the heavy quark mass expansion 
of the relativistic Bethe-Salpeter equation in [5] predicted a lower value 2.369 GeV of + 
mass which is still 50 MeV higher than the experimental data. Bardeen, Eichten and 
Hill interpreted them as the cs (0 + , 1 + ) spin doublet as the parity conjugate states of the 
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(0~, 1 _ ) doublet in the framework of chiral symmetry^] 1 (See also Ref. [7]). A quark- 
antiquark picture was also advocated by Colangelo and Fazio [9], Calm and Jackson [10], 
Godfrey [11]. Based on such a "conventional" picture, the various decay modes were 
discussed in Refs. [6, 9, 11]. 

Apart from the quark-antiquark interpretation, D s j(2317) was suggested to be a four- 
quark state by Cheng and Hou [12] and Barnes, Close and Lipkin [13]. Szczepaniak even 
argued that it could be a strong Dir atom [14]! But we think it would be very exceptional 
for a molecule or atom to have a binding energy as large as 40 MeV. 

Van Beveren and Rupp [15] argued from the experience with ao//o(980) that the low 
mass of _D s j(2317) could arise from the mixing between DK continuum and the lowest 
scalar nonet. In this way the + cs state is artificially pushed much lower than that 
expected from quark models. 

The recent lattice calculation suggests a value around 2.57 GeV for + state mass 
[16], much larger than the experimentally observed _D s j (231 7) and compatible with quark 
model predictions. The conclusion in Ref. [16] is that -D s j(2317) might receive a large 
component of DK and the physics might resemble ao//o(980). Such a large DK compo- 
nent makes lattice simulation very difficult. 

In this paper we shall use QCD sum rules [20] in the framework of the heavy quark 
effective theory (HQET) [21] to extract the masses since HQET provides a systematic 
method to compute the properties of heavy hadrons containing a single heavy quark via 
the 1/rriQ expansion, where ijiq is the heavy quark mass. The masses of ground state heavy 
mesons have been studied with QCD sum rules in HQET in [22, 23, 24]. In [17, 18, 19] 
masses of the lowest excited non-strange heavy meson doublets (0 + , 1 + ) and (1 + , 2 + ) were 
studied with the sum rules in HQET up to the order of O{l/mo). These masses were 
also analyzed in the earlier works [25] with sum rules in full QCD . In this short note we 
extend the same formalism in [17, 18] to include the light quark mass in calculating the 
cs mesons. We shall also use more stringent criteria for the stability windows of the sum 
rules in the numerical analysis. 



2 rriM — rriQ at the leading order of HQET 

The proper interpolating current J"pj" J for a state with the quantum numbers j, P, je in 
HQET was given in [17]. These currents were proved to satisfy the following conditions 

<<> ./-Krai)/. />'.,;) = f Pjl 5 jr 5 PP ,5 jd ,f^ , (i) 

i (o|t (./^-•i.r)./:/,;: •((.)) |o> = m^W- 1 ^^ • • 

x J dtS(x - vt) U PJe (x) (2) 

in the itlq — > oo limit, where r] ai " aj is the polarization tensor for the spin j state, v 
is the velocity of the heavy quark, g^ = g af3 — v a v 13 is the transverse metric tensor, S 

1 The existence of parity doublets has also been shown by combining chiral symmetry and heavy quark 
symmetry in Bethc-Salpeter approach in [8]. 
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denotes symmetrizing the indices and subtracting the trace terms separately in the sets 
(«i ■ • -aj) and ■ • • /3j), fpj e and Hpj e are a constant and a function of x respectively 
which depend only on P and jg. 
We consider the correlator 

nH = i J d 4 xe^(0\T (x)j]^r (0)) |0) (3) 

where u> = 2v ■ k. It can be written as 

f 2 

n(w) = - ' Je + higher states , (4) 

2A j,PJe - u 

Where Aj,p,j l = limm^^mMj P H — ttiq)- On the other hand, it will be calculated in 
terms of quarks and gluons. Invoking Borel transformation to Eq. (3) we get 

fp • e t = — / p(uj)e~T + condensates . (5) 

,M 7T J2m q 

where m q is the light quark mass, p{u) is the perturbative spectral density and u c is 
the threshold parameter used to subtract the higher state contribution with the help of 
quark-hadron duality assumption. 

We shall study the low-lying (0 + , 1 + ), (l + ,2 + ) (cs) states. The values of various QCD 
condensates are 

(ss) = -(0.8 ± 0.1) * (0.24 GeV) 3 , 
(a s GG) = 0.038 GeV 4 , 

m 2 = 0.8 GeV 2 . (6) 

We use m s (lGeV) = 0.15 GeV for the strange quark mass in the MS scheme. We use 
Aqcd = 375 MeV for three active flavors and Aq C d = 220 MeV for four active flavors. 
The sum rules with massless light quarks have been obtained in [18, 17]. 

In the numerical analysis of the QCD sum rules we require that the high-order power 
corrections be less than 30% of the perturbation term. This condition yields the minimum 
value T min of the allowed Borel parameter. We also require that the pole term, which is 
equal to the sum of the cut-off perturbative term and the condensation terms, is larger 
than 60% of the perturbative term, which leads to the maximum value T max of the allowed 
T. Thus we have the working interval T min < T < T max for a fixed u c . If T min > T max , 
we are unable to extract useful information from such a sum rule. In the ideal case, 
the difference between the meson mass mu and the heavy quark mass mg (or other 
observables) is almost independent on T for certain values of uo c . Namely, the dependence 
on both the Borel parameter and continuum threshold is minimum. In realistic cases, the 
variation of a sum rule with both T and u c will contribute to the errors of the extracted 
value, together with the truncation of the operator product expansion and the uncertainty 
of vacuum condensate values. 

In subsections 2.1 and 2.2 we give sum rules for (0 + , 1 + ) and (l + ,2 + ) doublets using 
interpolating currents with derivatives, respectively. For comparison we discuss sum rules 
for (0 + , 1 + ) using interpolating currents without derivatives in subsection 2.3. 
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Figure 1: The variation of A(| ) (in unit of GeV) of the (0 + , 1 + ) doublet with T and 
uj c for the derivative currents. The vertical and horizontal axes correspond to A and T. 
From top to bottom, the curves correspond to uj c being 3.1, 2.9, 2.7 GeV respectively. 



2.1 (0 + , 1 + ) Doublet With Derivative Currents 

We employ the following interpolating currents [17] 



1 v 

'0.+.2 " ^= 

j[% = -L h vl 5 ^(-i) p t s . (7) 



1 T 



where = D^ 1 — (v ■ D) v^, with = — ig A^ 1 being the gauge- covariant derivative. 
We obtain the sum rule relevant to A 

f'l i/ 2 e" 2A/T = r + 2™^ 3 - 6mV - 12m 3 s u)e-^ T du 

' 1 6A7T 2 J2m s 

~ < + \rril {as) - ^ (a s G 2 ) . (8) 

Taking the derivative of the logarithm of the above equation with respect to (™) one 
obtains the sum rule for A. Substituting the obtained value of A in (8) one obtains the sum 
rule for f' + jy 2 . In Figure 1, the variation of A(l/2 + ) with T and uo c is shown. According 
to the criteria stated above the working range is 0.38 < T < 0.58 GeV. We have 

A(l/2 + ) = (0.86 ±0.10) GeV, 

f' +>1/2 = (0.31 ±0.05) GeV 5/2 , (9) 
where the central value corresponds to T = 0.52 GeV and u c = 2.9 GeV . 
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2.2 (l + ,2+) Doublet 

For the (1 + , 2 + ) doublet, by using the following interpolating currents [17] 

4,+;r = /| (V a Q2 + tp - § ^r 2 a) <? , (io) 

we obtain the sum rule 

/+, 3 / 2 e- 2A/T = -^~ 2 f C + 2m s uj 3 c - 6m^ 2 - 12m 3 a u]e^ T 

mg (ss) - h^G*)T + -L^s) - ^r(^G 2 ) . (11) 

^From the sum rule, the variation of A(3/2 + ) with T and uj c is plotted in Figure 2, . 
The working range is 0.55 < T < 0.65 GeV. We have 

A(3/2 + ) = (0.83 ±0.10) GeV, 

/ +l3 /2 = (0.19 ±0.03) GeV 5/2 , (12) 
where the central value corresponds to T = 0.62 GeV and u c = 3.0 GeV . 
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Figure 3: The variation of the ratio of the pole contribution and the perturbative piece of 
the (0 + , 1 + ) doublet with T and uj c for the non-derivative currents. From top to bottom, 
ijj c = 3.1, 2.9, 2.7 GeV respectively. 



2.3 (0+, 1+) Doublet With Currents Without Derivative 

Possible different currents for (0 + , 1 + ) doublet are the non-derivative currents 

Jl,+,2 = -j= Krii s ■ 

With the non-derivative currents the sum rule reads 

3 



(13) 



,2 -2A/T _ 

/+,l/2 e " 167r 2 



[u 2 - 2m s co - 2m 2 s }e~ uj/T du 
1 



2m, 

1 ,_ v rn s ._ . 
+ 2^ + 4T (s_s) 8T* 



m 2 (ss) . 



(14) 



Requiring that the condensate contribution is less than 30% of the perturbative term, 
we get T m j n = 0.75 GeV. In Figure 3, we show the ratio of the pole term, i.e the sum of 
the cut-off perturbative term and condensation terms, for the sum rule (14). In the whole 
range of T > T m j n , this pole contribution is less than 40%. Hence, There is no stability 
window for this sum rule satisfying our criteria stated before. 

If we arbitrarily loosen the analysis criteria and require the pole contribution is greater 
than 30% only, we get the working range 0.75 < T < 1.2 GeV. In Figure 4, the variation 
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of A(l/2 + ) with T and uo c is shown. Numerically, 

A(l/2 + ) = (1.30 ±0.15) GeV, 

/ +)1/2 = (0.39 ± 0.05) GeV 3 / 2 , (15) 

where the central value corresponds to T = 1.0 GeV and u> c = 2.9 GeV. Because of 
the weaker criteria used these results are less reliable and will not be used in the final 
numerical results. 

Here we would like to make some remarks. Usually, the currents with the least number 
of derivatives are used in QCD sum rule approaches. The sum rules then are less sensitive 
to the threshold energy uj c . However, it is pointed out in [17] that in the non-relativistic 
limit the coupling constant of these currents to the P wave states vanishes. If this coupling 
constant is suppressed due to this reason, the relative importance of the contribution of 
the DK and other states in continuum in the sum rules which are usually neglected 
would be enhanced. Besides, it is shown in [26, 27] by using soft pion theorem that the 
contribution of the Dir continuum is large in the sum rule with the non-derivative current 
for the + state of the non-strange D system and significantly decreases the value of A. 
Similar method of calculation is not good in the present case, but it indicates that the 
contribution of the DK continuum with the non-derivative current may be large too. 
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3 Sum rules at the 1/mg order 

To the order of 1/mg, the Lagrangian of HQET is 

£ cS = h v iv-Dh v + — JC + —— 5 + 0(l/mn) , 
2m Q 2m Q w 



(16) 



where h v (x) is the velocity-dependent effective field related to the original heavy-quark 
field Q(x) by 



h v (x) = e m ^ x ^^Q(x) 



(17) 



/C is the kinetic operator defined as 



K, = h v {iD t ) 2 h v , 



S is the chromo-magnetic operator 
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<S = -C mag (m Q / '//) h v a^G^h v , 



(19) 



where C mag = (^j^j ^ A) = 11 ~ 2n,/3. 



Considering C(l/mg) corrections the pole term of the correlator on the hadron side 
becomes 



U{u) 



pole 



(f + Sff 



25mf 2 + 2f5f 



2(A + <5m) - 2A - (2A - cu) 2 2A - oo ' 



(20) 



where 5m and 5f are of order O (1/mg). 

To extract 5m in (20) we follow the approach of [23] to consider the three-point corre- 
lation functions 



'\M h r:, (u, oj ) = / 



2 j d 4 xd*ye ik -*- ik '-v(0\T(j^p(x) O(0) J^(y)) |0> , (21) 



where O = JC or S. The scalar function corresponding to (21) can be represented as the 
double dispersion integral 

1 f p (s, s')dsds' 



6011(00,00') = — J 



(s - u)(s' - u') 



The pole parts of 5oTl(u, u>'),0 — 1C,S are 



5 K \l(uj,uj') pole 



5 s Il(uj,uj') po i e 



f 2 G K (iu>) + f 2 G K (uj) 



(2A-w)(2A-w') 2A — uj 2k- u' ' 



d M f 2 Z 2 
77 + «mJ 



(2A-u;)(2A-u;') 



2A - 2A - u' 



(22) 



(23) 
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where 

K hP,h = (j, p Ji\h v (iD ± ) 2 h v \j,PJi) , 

2d M T ljiP:jl = (j, P, ji\h v gcr^G^hvlj, P, ji) , 

dM = d j:jl , = 2ji + 2, d jl+l/2ijl = -2j t . 

Letting u = uj' in eqs. (22) and (23) and comparing with (20), one obtains [23] 

The single pole terms in (22) and (23) come from the region in which s(s') = 2A and 
s'(s) is at the pole for a radial excited state or in the continuum. They are suppressed 
by making the double Borel transformation for both uj and uj'. The Borel parameters 
corresponding to uj and uj' are taken to be equal. One obtains thus the sum rules for K 
and S as 

f Ke -2A/r = r u ° r° duJ djuj'e-^/ 2T p K {uj,J) , (25) 

J 2m s J 2m s 

/2 Se -2A/T = r r fadJe-WVVTp^J) . (26) 

J2m s J2m s 

In this section we shall neglect the m s corrections to K and S. We obtain for the 
1 + 

jf = - doublet with derivative currents 
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/; 2 1/2 E e- 2A/T = f* f »') e-^dsds' + -i-(a.GG) T 3 , (27) 

' J2m 3 J2m s 48/T 



with 



-°( s > s ') = ^^^^^'{^(^ " - s') + (s <-> s')} . 

yo 7T° 



We can eliminate the dependence of the K and E on /+,i/2 and A through dividing the 
above sum rules (27) by the sum rule (8). We use the same working windows as those of 
two-point sum rules for A in evaluating 0(1/ mq) corrections. Numerically we have 

K 1/2 = (-1.60 ±0.30) GeV 2 

Iq /2 = (0.28 ±0.05) GeV 2 . (28) 

1 + 

The variations of K and £ with T and u) c for the jf = - doublet are shown in Figure 5 

2 

and 6, respectively. 

3 + 

For jf = - doublet, we have 

f 2 w3/2 Ze- 2i ' T = f p(s,s') e-'&dsds' + -i-(a.GG) T 3 . (29) 





Figure 6: The variation of £1/2 (in unit of GeV 2 ) of the (0 + , 1 + ) doublet with T and uj c 
for the derivative currents. From top to bottom, uj c = 3.1, 2.9, 2.7 GeV respectively. 
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Figure 7: The variation of K^/ 2 (in unit of GeV 2 ) of the (1 + , 2 + ) doublet with T and u c . 
From top to bottom, u c = 3.2, 3.0, 2.8 GeV respectively. 



with 

Pi*, o = wP^ Cmaa{s ' A{s ~ r Ms ~ s ' ] + (s ~ s ' )] ■ 

3 + 

The variations of K and S with T and u; c for the jf = — doublet are plotted in Figure 
7 and 8 respectively. Numerically we have 

K 3/2 = (-1.64 ±0.40) GeV 2 

S 3/2 = (0.058 ± 0.01) GeV 2 . (30) 

The spin-symmetry violating term S not only causes splitting of masses within the 
same doublet, but also causes mixing of states with the same j, P but different ji. In Ref. 
[18] corrections from the mixing is found to be negligible. We omit this effect here. 

4 Numerical Results and Discussions 

We present our results for the cs system assuming the HQET is good enough for excited 
D mesons. For the doublet (0 + , 1 + ) with the derivative current, the weighted average 
mass is 

j {m D * o + 3m D *J =m c + (0.86 ± 0.10) + — [(0.40 ± 0.08) GeV 2 ]. 
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Figure 8: The variation of £ 3 / 2 (in unit of GeV 2 ) of the (1 + , 2 + ) doublet with T and u c . 
From top to bottom, u c = 3.2, 3.0, 2.8 GeV respectively. 



The mass splitting is 

mn tl - m D * = ±- [(0.28 ± 0.05) GeV 2 ] . 
For the (1 + , 2 + ) doublet we have 

7^ (3m Dsl + 5m D *J =m c + (0.83 ± 0.10) + — [(0.41 ± 0.10) GeV 2 ] . 

(31) 

The 1 + , 2 + mass splitting is 

m D * -m Dl = — [(0.116 ± 0.06) GeV 2 ] . 

32 m c 

The results for bs system are obtained by replacing m c by m& and multiplying S by 
0.8 (since C mag ~ 0.8 for B system by using the values of Aq CD given in Sect. 2) in above 
equations. 

Choosing m c to fit the experimental value 

~ (3m Dsl + 5m D . ) = 2.56 GeV, 

where we again neglect the mixing between two 1 + states, we obtain m c = 1.44 GeV. 
Using this m c value we obtain the following numerical results. The 1 + , 2 + mass splitting 
is 

m D * 2 - m Dsl = (0.080 ± 0.042) GeV 
12 



which is consistent with the experimental value 37 MeV within the large theoretical 
uncertainty. Experimentally, the mass splitting in the (l + ,2 + ) doublet in D s system 
is almost equal to that in D system. This justifies neglecting the m s correction to the E 
term in our calculation. For (0 + , 1 + ) doublet, the weighted average mass is 

\ (m D * +3m D » ) = (2.57 ±0.12) GeV. 

A x sO sl / v f 

The mass splitting is 

m D * si - m D - to = (0.19 ± 0.04) GeV . (32) 

Therefore, the + mass is predicted to be m + = (2.42 ± 0.13) GeV. This is consistent 
with the experimental value 2.317 GeV, though the central value is 100 MeV larger than 
data. If the 1 + D s state of mass 2.460 GeV found by CLEO [2, 1] is assigned as the other 
member of the (0 + , 1 + ) doublet, then the observed mass splitting 0.143 GeV is consistent 
with (32). Notice that the £1/2 value in Fig. 6 changes very slowly between T = 0.35 — 0.8 
GeV. Therefore the result (32) is insensitive to the stability window used. 

We would like to note that the stability of the sum rules (in particular, those for the 
K and £) obtained by us is not as good as those for the ground states, as can be seen 
from the figures. And for the charm flavor hadrons, the l/m c corrections are significant. 
So the predictions on masses have a relatively large uncertainties, which are estimated by 
given errors. 

If we replace the strange quark condensate by up/down quark condensate and m s by 
the zero up/down quark mass, we can extract the non-strange excited D meson masses. 
For (l+,2+) doublet, the sum rule window is 0.57 < T < 0.67 GeV for cu c = 2.8 - 3.2 
GeV. The results are 

A 3/2 = (0.73±0.08)GeV, 
K 3/2 = -(1.6 ±0.4) GeV 2 , 

S 3/2 = (0.057 ± 0.01)GeV 2 (33) 

with the central value at T = 0.62 GeV. For (0 + , 1 + ) doublet, the sum rule window is 
0.46 < T < 0.6 GeV for oj c = 2.7 - 3.1 GeV. The results are 

A 1/2 = (0.79±0.08)GeV, 
K 1/2 = -(1.57 ±0.4) GeV 2 , 
E 1/2 = (0.286 ± 0.05)GeV 2 (34) 

with the central value at T = 0.53 GeV. The effect of strange quark is to make the mass 
of D s mesons a little bit larger, as expected. Experimentally, Belle collaboration found 
m Do = (2290 ± 22 ± 20) MeV and m D » = (2400 ± 30 ± 20) MeV recently [29] while 
CLEO collaboration found mc; = 2400^2 MeV [30]. The j 3 / 2 non-strange doublet mass 
is known precisely [28], m D2 = 2460 MeV and m Dl = 2420 MeV. Our results (33)-(34) 
are consistent with the experimental data within theoretical uncertainty. 

In summary, we have calculated the masses of the excited (0 + , 1 + ) and (l + ,2 + ) dou- 
blets for the cs system to the 1/mg order in the HQET sum rules. The numerical results 
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imply that the D sJ (2317) and £> sJ (2460) observed by BABAR and CLEO can be consis- 
tently identified as the (0 + , 1 + ) doublet with j t = | + within the theoretical uncertainty. 
Especially, the mass splitting in the (0 + , 1 + ) doublet is reproduced quite well. The re- 
pulsion between these states and the DK and DK* continuum may help to lower their 
masses. In the framework of the sum rule this effect should come from the contribution 
from the DK and DK* continuum to the dispersion integral. 
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SRF for ROCS, SEM. 
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